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We study braneworld models in the presence of scalar field in a five-dimensional geometry with 
a single extra dimension of infinite extent, with gravity modified to include a function of the Ricci 
scalar. We develop a procedure that allows to obtain analytical solution for the braneworld con¬ 
figuration in a diversity of models, in the much harder case where the Ricci scalar is non constant 
quantity. 

PACS numbers: 11.27.+d 


I. INTRODUCTION 

The idea that our four-dimensional Universe may be 
viewed as a subspace in a higher-dimensional scenario has 
attracted the attention of the scientific community work¬ 
ing in different fields, from elementary particle physics 
to cosmology, over the last fifteen years. In the first 
consistent models, space-time was described as a five¬ 
dimensional anti de Sitter ( AdS 5) geometry, with a sin¬ 
gle extra spatial dimension of infinite extent pQ. This 
is known as the thin braneworld scenario. This setting 
was soon modified to give rise to a braneworld config¬ 
uration sourced by a scalar field, giving rise to what is 
generally known as the thick brane scenario EMU- The 
thick brane scenario has been explored from many differ¬ 
ent perspectives In particular, given our ignorance 

about the gravitational dynamics in higher-dimensional 
backgrounds, models of the F(R) type, where R is the 
Ricci scalar, have also been studied in several works, e.g., 
png. Some of these studies allowed us to test the ro¬ 
bustness of the thick braneworld scenario under mod¬ 
ifications of the gravitational dynamics away from the 
standard Hilbert-Einstein action. Exploration of linear 
stability against tensorial perturbations is then possible 
if we have control on the background solutions. It is for 
this reason that finding exact solutions of different con¬ 
figurations becomes a useful tool to test the robustness 
of the thick braneworld scenario. 

In this work, we follow the lines of m and consider 
braneworld models described by theories of gravity of 
the F(R) type, but we focus on the general case of non¬ 
constant scalar curvature R. For the sake of generality, 
we take the scalar field with standard and generalized 
kinematics, but we start with the simpler possibility, with 
scalar field models having standard kinematics. Despite 
the freedom allowed by this approach, we show that it 
is possible to find analytical solutions for several specific 
models using a well defined procedure which consists on 
requiring that the warp factor be proportional to a hy¬ 
perbolic function of the extra dimension. The possibility 
to obtain analytical solutions with this approach allows 


us to study the stability of the gravity sector. In particu¬ 
lar, we find that the models considered are robust against 
fluctuations in the metric. 

To make the investigation easier to understand, we 
organize the subject of the work as follows. In Sec. II 
we investigate the braneworld model with a polynomial 
modification of the Einstein-Hilbert action on general 
grounds, with the source scalar field having standard 
kinematics. We then study several specific models to 
illustrate that the procedure to generate analytical so¬ 
lutions works nicely. In Sec. Ill we investigate the sta¬ 
bility of the models worked out in Sec. II, and verify 
their robustness against fluctuations in the metric sec¬ 
tor. In Sec. IV we extend our methodology to models in 
which the scalar field has generalized kinematics. This is 
much harder, but we show that the procedure still works, 
leading to the construction of exact solutions for several 
specific models. We end the work in Sec. V with some 
comments and conclusions. 

II. BRANEWORLD MODEL 

We start with a 5 D model, with F(R) gravity coupled 
to the real scalar field (f> in the form 

s = jd 5 xV\i\ + • (i) 

We take 47 tG^ 5 ) = 1 and use g = det(g a b). We investigate 
the case of a flat brane, with the line element given by 
ds 2 = e 2A ?7 Ml/ — dy 2 where e 2A is the warp factor, 77^ is 
the 4-dimensional Minkowski metric, and y = x 4 is the 
extra dimension. 

We suppose that both A and 4> are static and depend 
only on the extra dimension, such that A = A(y) and 
(f> = <j>(y). In this case, the variation of the action 0 
over the scalar field gives the following equation 

<tt' + 4 A' 4 ! = , (2) 

where prime denotes derivative with respect to the extra 
dimension, and = dV/dcj). In this case we obtain the 


2 


scalar curvature as 


R = 20A ' 2 + 8 A" . 


( 3 ) 


The equations obtained by the variation over the metric 
can be combined to obtain the following equations 

<\J 2 = - | A"F r +4 (5 A' 2 A"- 5 A" 2 - A A'A’"- A"") F RR - 
-32 (5 A! A" + A'") 2 F rrr , (4a) 

V{4>) = -\F{R) + i (5A" + 8 A 12 ) F r - 

-2 (35 A 12 A" + 12 A'A" 1 + 5A" 2 + A"") - 

-16(5 AA"+ A!") 2 F rrr . (4b) 


Also, the energy density of the scalar field can be written 
in the general form 


p(y) = e 2A 


-±F(R) + ^(A" +4A i 2 )F r - 

-4 (15 A ' 2 A"+ 5 A" 2 + 8 A' A"'+ A"") F RR - 


-32 (5A'A" + A'") 2 F rrr 


( 5 ) 


For concreteness, in this work we choose to work with 
the family of gravitational Lagrangians 

F(R) =R + aR n , (6) 

where a is a real number and n is a positive integer. With 
this choice, the equations Q become 

4 (/ 2 = -^A" - ^A" (20A ' 2 + 8 A ")”' 1 + 

+4an(n— 1) (5A ,2 A"- 5A" 2 - 4A'A'"-A"") x 
x (20A' 2 + 8A") n ~ 2 — 32 an(n - 1 ){n - 2) x 
x (5A'A" + A'") 2 (20A' 2 + 8A")” -3 , (7a) 

V{<j>) = -3A /2 - ^A"— ^a ( 20 A ,2 + 8 A") ti - 

+ ^ (5A" + 8 A' 2 ) (20A ' 2 + 8 A ")” -1 - 
-2cm(n-l) (35A , 2 A" + 12A , A" , +5A" 2 +A"") x 
x (20A 72 + 8 A") n 2 — 16an(n— 1)(n—2) x 
x (5A , A"+A "') 2 (20A 72 + 8 A ")” -3 . (7b) 


The energy density then turns into 


p(y) = e 2A [ - 3A ' 2 - ^A" - ^a (20A 72 + 8 A") n + 


+ °^ (A" + 4A' 2 ) (20A ' 2 + 8 A") n_1 - 
-4an(n- 1) (15A 72 A"+5A" 2 +8A 7 A" 7 +A 7777 ) x 
x (20A ' 2 + 8 A") n ~ 2 - 32 an{n - l)(n - 2) x 


x (5A'A" + A "') 2 (20A ' 2 + 8 A ")” -3 


( 8 ) 


The challenge now is to find a consistent solution of the 
gravitational and scalar-field equations given above. To 
proceed, we assume that the space-time behaves asymp¬ 
totically (when y —> oo) as an anti-de Sitter (AdS) geom¬ 
etry. This requirement can be implemented by choosing 
the warp function in the form 

A{y) = B In (sech(fcy)) , (9) 

where B and k are real and positive parameters. The 
positivity of the parameters is necessary for the graviton 
to be localized on the brane in the thin brane limit |3 El- 
We will next show that the ansatz ^ is good enough to 
allow us to obtain analytical solutions in several examples 
of current interest. 


A. First Example 


As a first example we consider a quadratic term in 
equation © as done in [19] . We take F(R) = R + aR 2 , 
and write the equations ([7]) as 

<t> r2 = - 1 A" - 4a (l 6 A" 2 + 5 A' 2 A"+ 8 A'A'"+ 2 A ""), (10a) 

y (</>) = — ^ A" — 3 A ' 2 (10b) 

-2a (8A" 2 +69A' 2 A"+ 10A' 4 +24AA"'+2A""). 
Also, the energy density ([ 8 ]) is simplified to the form 


p(y) = -e 2A 


3A ' 2 + ^A" + 


( 11 ) 


+4a (5A ' 4 + 37 A' 2 A" + 12A " 2 + 16A'A" 7 +2 A"") 


Since the scalar field is real, the parameters k and B 
characterizing the warp factor and the scale a in the La- 
grangian must satisfy the following constraint 


-3 


8 fc 2 (8 + 16B + 5B 2 ) 


< a < 


32fc 2 (l + 4B)' 


( 12 ) 


We can now substitute equation © in (jTOaJ) to obtain 
the differential equation for the field (f> as 


<j)' 2 = a\sech 2 (ky) — a 2 sech 4 (fcy), 


(13) 


where ai and a 2 are constants, given by 

ai = ^Bk 2 +4aBk 4 (8 + 16B+ 5B 2 ), (14a) 

a 2 = 4aBk 4 (6 + B)(2 + 5B). (14b) 


The issue now is to obtain analytical solution of the equa¬ 
tion (13) with a suitable choice of the parameters. For 


example, if we choose a = —3/[8fc 2 (8 + 16 B + 5 B 2 )] we 
get from (13) the new equation 


k I 6B(6 + B)(2^B) 

P 2 V 8 + 16B + 5B 2 ( V) ’ 


(15) 
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FIG. 1: Energy density ( |18[ ) depicted for k = 1, B = 1 and 
a = 0 (dashed line) , a = a s /2 (thicker line), a = a s (thick 
line) and a = 3a a /2 (thinner line), with a s = 3/250. 


which can be solved to give 


B. Second Example 

Let us now analyze the model with cubic term in the 
Ricci scalar, that is, F(R) = R + aR 3 . In this case, the 
above equations 0 become 

<// 2 = -1 A"+ 24a (25A ,4 A"~ 320A ,2 A" 2 - 80 A' 3 A'"- 
-52A" 3 - 112A'A"A'" - 8 A'" 2 - 
-20 A' 2 A'"' - 8A"A"") , (21a) 

V{4>) = - A"-3A ,2 +4a (lOOA' 6 - 1845A ,4 A"-92A" 3 - 
— 1584A ,2 A" 2 - 720 A' 3 A 1 " - 528 A 1 A"A 1 ” - 
—24 A'" 2 - 60 A' 2 A"" - 24A" A"") , (21b) 

and the energy density can be written as 


</>(y) = ±2 


1 6 B(6 + B)(2 + 5B) 


8 + 16 B + 5 B 2 


Back to equation (13), we note that the standard solution 


is obtained with a = 0; it is given by 


tanh (ky). (16) p{y) = -e 2A (s A' 2 + ^ A" - 8a(b0A' 6 -885A' 4 A”- 

-1272A! 2 A!’ 2 - 124 A" 3 — 480A' 3 A'"- 432 A'A"A'"~ 
—24A"' 2 - 60A' 2 A"" - 24A" A"") ) . (22) 


<j>{y) = \l arcsin tanh (ky) 


(17) 


Additionally, with the choice 0 can write the energy 
density 


as 


p(y) = —B 2 k 2 (3 + 20 aB 2 k 2 ) sech 2 B (ky) 
+aq (1+2 B) sech 2 B+2 {ky) — 

—02 (1+1?) sech 2B+4 (fcy), 


(18) 


where aq and 02 are given by (141. An important feature 


of this model is the presence of the split of energy density. 
This aspect was first pointed out in m , and it appears 
for 


a —> a x = 


9B + 3 

8 k 2 (49 B 2 + 60 B + 16) ' 


(19) 


This and (l2| ) impose that —5.68 < B < —2.56 or B > 
—0.25. In Fig. [l] we depict the behavior of the energy 
density for several values of a, using k = 1, B = 1 and 
a s = 3/250. 

Though for generic values of the parameters it is not 
always possible to get the scalar potential as a function 
of the field <j> explicitly, for the particular solution (16) it 
is indeed possible. In this case, 


VW = 


3k 2 B f 5B 2 + 24B + 12 


— fc 2 (l + 2B)4> 2 + 


4 \5B 2 + 16 B + 8 

k 2 / (1 + 2B)(8 + 16 B + 5 B 2 ) \ M4 
Y V B(6 + B)(2 + 5B) 


This potential engenders spontaneous symmetry break¬ 
ing and supports the kink-like structure already shown 

in (fT6l) . 


As we can see, the equation for the scalar field de¬ 
pends on the fourth derivative of the warp function, and 
suggests that we can use the above procedure to obtain 
exact solutions. We then suppose that the warp function 
is again given by §. With this, the differential equation 
for the scalar field can be written as 

cj)' 2 = bisech 2 (ky) — b 2 sech 4 (ky) + b 3 sech 6 (ky) , (23) 

where the constants b±, 62 and 63 are given by 
3Bk 2 

h = —— + 120aR 3 fc 6 (16 + 32 B - 5 B 2 ) , (24a) 

62 = 48aB 2 fc 6 (2 + 5£?)(16 + 66B — 5B 2 ) , (24b) 

b 3 = 24aB 2 k 6 (20-B)(2 + 5B) 2 . (24c) 

In this new case, there exist several possibilities to find 
analytical solutions, including © when a = 0 , as ex¬ 
pected. Taking B = 20 we see that b 3 = 0, which pro¬ 
vides the first non-trivial simplification. Assuming now 
that a is such that b\ = 0 (ie, a = l/(43008000fc 4 )) we 
can write the differential equation (23) as 


H = ±k 


4233 , 2 

- sprh 

140 


sech (ky). 


(25) 


This differential equation admits a solution of the form 

(26) 


<j>(y) = ±J——tanh(ky ), 


140 


(20) which has the same behavior as the solution (16) obtained 


in the previous model. In this case, the potential as a 
function of </> is given by 

V^) = ————— k 2 - 41k 2 Y + —k 2 ^. (27) 

yr> 840 v 4233 v ; 


840 
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FIG. 2: The solution (261 (solid line) and (291 (dashed line), 
depicted for k = 1. 



FIG. 3: The energy density ( |30[ ) depicted for k = 1, B = 20 
and a = 0 (dashed line), a = as a /2 (thicker line), a = a 3 S 
(thick line) and a = 3a3 S /2 (thinner line). 


Another possibility is obtained if 5 = (33 + -\/ll69)/5 
and a is such that b i = 0. In this case the differential 
equation (23) becomes 


4>' = sech 3 (ky) , (28) 


which leads to the new solution 


4>{y) = ±\j ^163+5-\/il69^ arctan^tanh(/cy/2) 


+ -sech(fcy) tanh(fcy) 


(29) 


In this case, we cannot get the scalar potential in terms of 
the field <j> explicitly. In Fig. [2] we illustrate the behavior 
of the solutions (26) (solid line) and (29) (dashed line), 
depicted for k = 1. In the general case, we can express 
the energy density as 


p(y) — ~B 2 k 2 (3 — 400 aB 4 k 4 ) sech 2B (ky) + 

+ 6 i(l + 25) sech 2B+2 {ky) + 

+6 2 (l+5) sech 2B+4 (ky) + 

+ ^(3+25) sech 2B+6 {ky). (30) 

Again, we can show that the energy density begins to 
split when 


95 + 3 

a -> a 3s = - 64Bfc4 + 3515 + 84 ) ' 


(31) 


C. Third Example 


Let us now study the case with 5(5) = R + aR 4 . Here 
the differential equation (7a) for the scalar field becomes 


= - 6 aA" (20 A' 2 + 8A") 3 + 

+48a (5 A' 2 A"- 5A " 2 -4 A'A'"-A"") x 
x (20A ' 2 + 8 A") 2 - 768a x 
x (5 A'A" + A"’) 2 (20A ' 2 + 8 A") . (32) 

Thus, using equation (|9|, we find that 

4>’ 2 = ci sech 2 {ky) + c 2 sech 4 (fcy) + 

+c 3 sech 6 (fey) + C4 sech 8 (fey), (33) 


where the constants are given by 
3 Bk 2 

ci = —+9600afc 8 5 5 (8 + 165-55 2 ) , (34a) 

c 2 = 1920afc 8 5 4 (2 + 55) (155 2 -1025 - 32) , (34b) 
c 3 = 1152afc 8 5 3 (2+55) 2 (8 + 525 - 55 2 ) , (34c) 
c 4 = 384a/c 8 5 3 (2+55) 3 (5 - 14). (34d) 


In this case we suppose that 

a —> a 1 = l/(64005 4 fc 6 (55 2 - 165 - 8 ), (35) 


then ci = 0. Also, we take 5 as one of the possible 
solutions of the following equation 


3 (8 + 525 - 55 2 ) 
2(5 - 14) 


= 


15B (15B 2 — 102B—32) 


(5 - 14) 


(36) 

With this restriction, the differential equation (33) be¬ 


comes 


We see that in this new situation the splitting appears for 
a negative. In Fig. [3] we show the behavior of the energy 
density near the splitting. Here we are using k = 1 and 
5 = 20, with a 3s = -183/187781120. 


*' = ± / -16B 560113 ikV) h s » h 2 W) • 

(37) 
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FIG. 4: The solutions (|39|), depicted for k = 1. 


where 

_ 3 (8 + 52B — 5B 2 ) 

° 2(2 + 5B)(B — 14) 1 j 

For equation ( [37] ) to be true, we must impose that c 4 > 0, 
with a = ai, which causes8/5 — 2v / 26/5 < B<— 2/5, or 
0< H <8/5+2v / 26/5, or B > 14. Solving the differential 
equation ( [37] ) we obtain, for the allowed values of B , the 
solution 

</>(y)=± (5.99 — 0.87 sech 2 (fct/))tanh(fcy), (39) 

for H = 17.448. It is depicted in Fig. [4] 

The energy density can be written as 

p(y) = -3B 2 k 2 (l - 8000 aB 6 k 6 ) sech 2 B (ky) + 
+ci(l+2S) sech 2 B+2 {ky) + 

+c 2 (1+B) sech 2B+i (ky) + 

+ +3+2 B) sech 2 B+6 (ky) + 

O 

+ ^(2+H) sech 2B+8 (fcy). (40) 

We can also see that in this case the energy density starts 
to split for a given by 


— 


1 + 3 B 


1024H 2 fc 6 (40 + 1765 + 1895 2 ) ' 


(41) 


To see how the splitting appears, in Fig. [5] we depict 
the energy density (40) for some values of a , with k = 1, 


B = 0.239 and a 4s = 0.00032. 

The above investigations show that the procedure is 
very efficient to obtain analytical solutions for n = 2,3, 
and 4, and this is important progress. However, we note 
that when n increases, the number of free parameters be¬ 
comes severely constrained. In particular, if n is greater 
then 4 this makes it harder to find analytical solutions. 
We can circumvent the problem relying on numerical in¬ 
vestigation, an issue that is out of the scope of the current 
work. 

To complete our study, we must determine if the modi¬ 
fications of gravity introduced in the above examples can 
contribute to destabilize the brane. This is done in the 
next section, where we investigate their linear stability 
in the gravity sector (tensorial modes). 



FIG. 5: Energy density (401 depicted for k = 1, B = 0.239 
and a = 0 (dashed line), a = a.As/2 (thicker line), a = a 4s 
(thick line) and a = 3 a 4 S /2 (thinner line). 


III. STABILITY 


The investigation of linear stability of the generalized 
braneworld models presented above can be carried out in 
different levels of complexity, depending on whether one 
is interested in the scalar field sector and/or in the gravi¬ 
tational side. Here, for simplicity, we focus on purely ten¬ 
sorial perturbations which, as is well known, are gauge 
independent. We can thus assume that the brane remains 
static (Gaussian normal coordinates) and write the per¬ 
turbed line element as (here x stands for a vector in the 
(3,1) space-time) 


ds 2 = e 2A ^ 


Vu,v + h^ v (x, y) dx^dx” - dy 2 . (42) 


Up to first-order in the tensorial perturbation the 
field equations yield 


^ h ; v + 2A'h^-\e- 2A n^h l , v + 


Fr+ 


„-2 A ; 


-2/lr 


~e 2A (d^d a h a „+d„d a h atJ ,-d^d v h 

1 9 M a !/ P( x, y) = 13 A'P'(x, y) + P"(x, y) - 
*□ (A) P(x, y) + (A'h + \h!) F' r + l -A'h' + 

+2V*f+20T- [4 A' 2 +A" (43) 
where h = ftT^ and the function P(x , y) is defined as 


P{x,y)=[h"+5A l h l +e~ 2A (d^d^h^-D^h^FRR.iU) 


As usual, we consider the perturbations to be transverse 
and traceless, i.e., we take 


= 0 and h = 0 


(45) 


which reduces Eq. (43) to the form 


(-d 2 -4A , d y +e- 2A nW-^dy ')V = 0. (46) 
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We now introduce a new coordinate z via the transfor¬ 
mation 


dz = e A ^dy 


(47) 


such that the line element is now conformally flat, and 
decompose the modes in the form 

M*. *) = e- ip - x e- 3A ^ l2 F~ 1/2 h^(z ). (48) 

By doing this, the 4-dinrensional components of obey 
the Klein-Gordon equation and the metric fluctuations of 
the brane satisfy the Schrodinger-like equation 


d?_ 
dz 2 


U{z 


hfj,u — 


nv ) 


(49) 


where 


1 / d(ln Fr) \ 2 1 d 2 F R 
4 V dz ) 2 Fr dz 2 


(50) 


The stability of the above modes can be established by 
noting that one can write 


dz 2 


+ U(z) = S*S, 


where 


o __d 3 , l d(hiF R ) 

dz 2 z 2 dz 


(51) 


(52) 


This decomposition puts forward that the operator —d 2 + 
U(z) is non-negative, which guarantees that the gravity 
sector is linearly stable. We note that the stability be¬ 
havior in the gravity sector only depends on the warp 
factor since R is a function of A(z), explicitly given by 
R = 4e~ 2A ( z ) (3A% + 2 A zz ). In addition, the zero mode 
can be written as 


Vo(z) = N 0 e 3 A ( z )/ 2 Ti /2 . 


(53) 


For the polynomial family of gravity Lagrangians de¬ 
fined in equation the potential of stability ( |50| can 
be written in the specific form 


, 9,n 3 . cm(n - 1) 

I/ «=4- 4 - + 2- 4 “-3SfT2x: 


/ i ?"" 1 \ 

v 1 + anR ” _1 / 


3 At + 8 AlA zz - A 2 - 2 A Z A ZZZ - A 


+ - 


any 


{n - 1 )R n ~ 1 (3A 3 Z - A Z A ZZ - A z 


w( 


1 + ani ?"- 1 V 3A 2 +2A Z 

2(„-2)- an(n - 1>X "~ 1 ' 


1 + anR’ 


n—1 


(54) 


whose zero mode is readily seen from (53) and reads 


Vo(z) = N 0 e 3A( - z ^ 2 \/l+anR n ~ 1 . 


(55) 


The next step is to make the change of variable (47) for 
the choice of warp function introduced in This can 
be done in a general way so that it allows us to obtain the 
following relationship between the z and y coordinates 


z = 


sinh (ky) 


lFl Cr 


1 1 -B 3 . , 2 

2 ’ 2’ _ Smh 


(fcy)), (56) 


where 2^1 is hypergeometric function. Note that this 
equation can not be inverted analytically for arbitrary 
B, though it is possible for some particular choices. For 
instance, if B = 1 we obtain y = arcsinh(fc^:)/fc. 

With all these elements, we are ready to investigate 
the linear stability of the polynomial models presented 
above. 


A. Case n = 2 

For the first example, involving a quadratic term in the 
Ricci scalar, the stability potential (54) can be written 
as 


U(z)~- -A 2 Z + - As* — (- 


8 a 


e 2 A(z) +8a(3A 2 z +2A zz ), 

3Af + 8A 2 A zz — A" zz — 2A Z A ZZZ — A zzz 

a a it 3A 3 — A Z A ZZ — A zzz \ 2 
—64a 


,e 2A W +8a(3A2 +2A ZZ ), 


(57) 


In this case, the parameter B can take any positive 
value and, for simplicity, we set it to B = 1. This choice 
has the clear advantage of allowing to perform analyti¬ 
cally the inversion of equation (56), so that the stability 
potential becomes 


U{z)=- 


3 k 2 2 — 5 k 2 z 2 

(1 + k 2 z 2 ) 2 ~ 

56 ak A (l — 6 k 2 z 2 ) 

(1 + k 2 z 2 ) 2 [l + k 2 z 2 — 8ak 2 (2 — 5 k 2 z 2 )] 
3136 a 2 k 8 z 2 

(1 +k 2 z 2 ) 2 [l+k 2 z 2 — 8ak 2 (2 — 5 k 2 z 2 )} 2 


(58) 


In addition, the zero mode is 


Vo(z) = 




(1 + k 2 z 2 ) 


3/4 


(59) 


The potential of stability (|58| and the zero mode (59) 
are depicted in Fig. [ 6 ] and Fig. [7J respectively. In both 
cases we used values of a around the splitting value. It 
should be noted that as the splitting of the energy den¬ 
sity occurs, the form of the potential at y = 0 changes 
dramatically from a deep minimum to a local maximum. 
As a result, the ground state becomes less localized, with 
the top of the wave function being lower and its width 
larger. 
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FIG. 6: Stability potential (581 depicted for k = 1 and and 
a = 0 (dashed line) , a = a s /2 (thicker line), a = a s (thick 
line) and a = 3a s /2 (thinner line) with a s = 3/250. 



FIG. 7: Zero mode ( |59| ) depicted for k = 1 and and a = 0 
(dashed line) , a = a B /2 (thicker line), a — a B (thick line) 
and a = 3a s /2 (thinner line) with a 3 = 3/250. 


B. Case n — 3 


Let us now analyze the model with cubic corrections in 
the Lagrangian. In this case the stability potential (54) 
becomes 


U (z) = ^A l + ^A ZZ - 96a 


3A 2 +2A* 


Le 4j4 ( z ) +48a(3^1? + 2 A zz ) 2 J 


3A z -{-8A z A zz — A zz — 2A Z A ZZZ — A z 


+192ae 4A(z) ( 


3 A z — A Z A ZZ — A zzz 
e^ A (A+48a(3A 2 z + 2A zz ) 2 


■ (60) 


Though we found two exact solutions for this model, the 
corresponding values of B do not allow for an analytical 
inversion of equation (56). Nonetheless, the inversion can 
be carried out numerically, which is done in Figure [8j 
using k = 1, B = 20 (solid line) and B = (33 +-\/ll69)/5 
(dashed line), with a = 1/43008000. For this model, the 
zero mode (55) can be written as 


7 io(z) = e 3A (A/ 2 ^ji + 48ae -4A ( 2 ) (sA 2 z + 2 A z ^ , (61) 

whose numerical representation appears in Fig[9j It is ap¬ 
parent from this plot that the localization of the ground 
state around the brane is tighter for smaller values of B 1 


with the potential minimum being deeper and the local 
maxima higher. 



FIG. 8: Stability potential (60 1 depicted numerically for k = 
1, a = 1/43008000, and B = 20 (solid line) and B = (33 + 
VH69)/5 (dashed line). 



FIG. 9: Numerical representation of the zero mode (611 for 
k = 1, a = 1/43008000, and B = 20 (solid line) and B = 
(33+ %/1169)/5 (dashed line). 


C. Case n = 4 


For n = 4, the stability potential ( |54[ ) reads 

m+ = 9 4 2 3 768a(3Al + 2 A zz ) 2 

fol 4 Az + 2 “ e 6A (A + 256a(3A2 + 2 A zz ) 3 

3 At + 8 A 2 t A zz — A 2 — 2 A Z A ZZZ — A zz - 


+768a(3A z +2A zz ) 

xU- 


(3 A z — A Z A ZZ — A zzz ) 2 


L e 6A(z) +2 56a(3A2+2A„) 3 J 
768a(3^2 + 2A zz f 


e 6 A (A + 256a(3A 2 + 2A ZZ ) 3 1 ' 


(62) 


Similarly as in the example of n = 3, the analytical in¬ 
version of equation (56) is not possible. Thus, once again 
we resort to a computer analysis to study the behavior 
of the stability potential, which appears in Fig [lOj For 
this model, the zero mode takes the form 

r) 0 (z) =e 3A ( z)/2 l+256ae~ 6A (A (3A 2 Z + 2A ZZ Y , (63) 









































where 



FIG. 10: Stability potential (621, depicted for k = 1 and 
5 = 17.448 and a = 1.36 x 10 _1 . 


and is plotted in Fig EH 

As informed in the beginning of the section, we have 
studied stability against tensorial perturbations. With 
this, we can then investigate several other related issues, 
among them the presence of resonances, the Newtonian 
potential, the entrapment of fermions and gauge fields, 
and other forms of fluctuations, in particular the addition 
of scalar perturbations. These problems have been stud¬ 
ied in several different contexts, for instance, in [31SHT41- 
[25] and in references therein. The point here is that such 
investigations are in general more involved and would 
usually rely on numerical study, leading us to issues that 
are out of the scope of the current work. Instead of fol¬ 
lowing such possibilities, we now broaden the scope of 
the work extending the scalar field sector, allowing that 
the source field engenders generalized dynamics, an issue 
that we deal analytically in the next section. 



FIG. 11: The zero mode (631, depicted for k = 1, B = 17.448 
and a=1.36 x 10~ 12 . 


IV. EXTENDING THE SCALAR SECTOR 

The methodology presented above to obtain exact solu¬ 
tions for brane configurations governed by different grav¬ 
itational Lagrangians can also be extended to scenarios 
in which, in addition, the source scalar field has non¬ 
standard dynamics. For the sake of generality, we can 
consider models of the form 

S = Jd 5 xV\g\ (- \f(R) + G{X) (64) 


X = ^V a </>V a <(> , (65) 

and G(X) is an arbitrary function of X. Modifications 
of this type in the scalar sector have been considered in 
the literature in diverse contexts, such as in inflationary 
models (as a way to produce non-Gaussianities) and in 
relation with dark energy models. Modifying the form of 
the scalar sector necessarily has an impact on its solutions 
as compared to the case of canonical scalar fields. How¬ 
ever, as tensorial perturbations do not get mixed with 
scalar modes, the stability analysis performed in the pre¬ 
vious section remains valid mum 

In this section, therefore, we focus on the search for 
analytic solutions in models with non-canonical (or gen¬ 
eralized) kinetic term, putting forward in this way the 
versatility of our approach in different scenarios with 
modified dynamics. Rather than providing an exhaus¬ 
tive list of models and their solutions, we just point out 
that successful results can be obtained for models of the 
form 


N M 

F(R) = J2 a n R n ; G(X) = b mX m , (66) 

n= 1 m= 1 

though we will only present results for the illustrative 
case 


F(R) =R- aR 2 ; G(X) = X- aX 2 , (67) 


where a and a are real constants, and a > 0 is necessary 
to ensure hyperbolicity [25] of the equation of motion of 
the scalar field. For other values of N and M the results 
are very similar. 

For the above model with quadratic corrections, the 
field equations become 

(f> 12 +a(t> ,A = —-A"+ (68a) 

+4a (2 A""+16 A ,,2 + 8A'A'" + 5 A' 2 A”), 
R(0)-“<(/ 4 = -||a"-3A ,2 + (68 b) 

+ 2a (8A" 2 +69A' 2 A"+10A ,4 +24A , A , "+2A""). 


Note that if a = 0 and 


» — a we get back the the 
equations (fl0|), as expected. On the other hand, if o/ 


0, the equations involve up to the fourth power of cf/, 
which makes it more difficult to find solutions. However, 
as suggested previously, we can use the standard warp 
function ansatz ([£]) to see how it performs here. With 
that choice, equation (68a) turns into 


<// 2 +a<// =-(3 —232a) sech (t/) + 196a sech (y), (69) 


where we used 5=1. This equation can be solved if we 
assume a suitable choice for the parameters a and a. For 
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example, if a = 784a/(3—232a) 2 , we obtain the following 
solution 


/ 3 

-| 

w — — 116a arcsin 

tanh(y) 


(70) 


For the solution to be real we also must impose that 
a <3/232. Note that the asymptotic value of the solution 
is changed to </— >±4> with cf>= (-^2^/4 )-^3— 2320 . With 
the solution (Tfol) we can use equation (|68b|) to obtain the 


scalar potential as 


V(<j>) = -3 + 20a + ^ (3 - 232a) cos 2 (| t) + 


+ 343a cos ( — ■= 


V2 &) ' 


where we have V (</>) = —3+20a. In Fig. 


(71) 

12] w e depict the 
( |7ip . In the two 

;he value a = 0 is represented by the continuous 
line and a = 0.005 by the dashed line. We chose the value 
a = 0.005 because it is inside the allowed interval. 


solution ([70]) and in Fig. [T3|the potential 
figures, t 



FIG. 12: The scalar field solution (701, 
(solid line) and a = 0.005 (dashed line). 


depicted for a = 0 



FIG. 13: The scalar potential [71 [ 
line) and a = 0.005 (dashed line). 


, depicted for a 


0 (solid 


One can also express the energy density as follows 

p{y) = -{ 3-20a) S' 2 + ^(3-232a) S' 4 + 392a S 6 . (72) 

Note that in this case the split happens for negative val¬ 
ues of a, which are not allowed on physical grounds. In 



FIG. 14: The energy density ( |72[ ) , depicted for a = 0 (solid 
line) and a = 0.005 (dashed line). 


Fig. [14] we depict the behavior of the energy density for 
a = 0 (solid line) and a = 0.005 (dashed line). 

As stated before, the inclusion of the term X 2 does 
not contribute to destabilize the brane, and the stability 
behavior is similar to that obtained for the first example 
in Sec. Ill A. 


V. COMMENTS AND CONCLUSIONS 

In this work we have studied thick braneworld mod¬ 
els with generalized gravitational dynamics of the form 
F = R + aR n , and have found analytical solutions for 
n = 2,3 and 4. As explicitly shown, as n increases, 
the search for analytic solutions becomes more and more 
complicated. However, we have been able to show that 
parameterizing the warp factor as the hyperbolic secant 
of the extra dimension very much simplifies the equations 
and allows to find analytic solutions in explicit form. This 
procedure also works when the kinetic term of the scalar 
field sourcing the thick brane scenario is changed, which 
further supports the strength of our method. 

It should be noted that the field equations that deter¬ 
mine the form of the scalar field and its potential gener- 
ically involve up to fourth-order derivatives of the warp 
function, while our prototypical ansatz only depends 
on two free parameters. This means that there is still 
freedom to further generalize the warp function ansatz 
to include two more parameters. Though this possibil¬ 
ity might have some interest from a mathematical point 
of view, we believe that our choice with two parameters 
provides enough freedom to study the relevant physics 
of these models, namely, to show that consistent exactly 
solvable models exist and that they are stable against 
tensorial perturbations. 

The results presented in this work motivate us to sug¬ 
gest that one uses the same methodology to find exact 
solutions in other scenarios of current interest, with grav¬ 
ity modified as in the Born-Infeld, Gauss-Bonnet and in 
other cases of current interest. 
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